Atmos. Chem. Phys., 13, 2423-2434, 2013 
www.atmos-chem-phys. net/1 3/2423/20 1 3/ 
doi: 10. 5 194/acp- 13-2423-20 13 
© Author/ s) 2013. CC Attribution 3.0 License. 



Atmospheric f 

Chemistry | (Rfcf 
and Physics 1 


Black carbon vertical profiles strongly affect its radiative 
forcing uncertainty 

B. H. Samset 1 , G. Myhre 1 , M. Schulz 2 , Y. Balkanski 3 , S. Bauer 4 , T. K. Berntsen 1 , H. Bian 5 , N. Bellouin 6 , T. Diehl 7 - 8 , 
R. C. Easter 9 , S. J. Ghan 9 , T. Iversen 2 - 11 - 16 , S. Kinne 10 , A. Kirkevag 2 , J.-F. Lamarque 12 , G. Lin 13 , X. Liu 9 , 

J. E. Penner 13 , 0. Seland 2 , R. B. Skeie 1 , P. Stier 14 , T. Takemura 15 , K. Tsigaridis 4 , and K. Zhang 9 10 

'Center for International Climate and Environmental Research - Oslo (CICERO), Oslo, Norway 
2 Norwegian Meteorological Institute, Oslo, Norway 

3 Laboratoire des Sciences du Climat et de l’Environnement, CEA-CNRS-UVSQ, Gif-sur- Yvette, France 

4 NASA Goddard Institute for Space Studies and Columbia Earth Institute, New York, NY, USA 

5 Joint Center for Earth Systems Technology, University of Maryland Baltimore County, MD, USA 

6 Met Office Hadley Centre, Exeter, UK 

7 NASA Goddard Space Flight Center, Greenbelt, MD, USA 

8 Universities Space Research Association, Columbia, MD, USA 

9 Pacific Northwest National Laboratory, Richland, WA, USA 

10 Max Planck Institute for Meteorology, Hamburg, Germany 

11 Department of Geosciences, University of Oslo, Oslo, USA 

12 NCAR Earth System Laboratory, National Center for Atmospheric Research, Boulder, CO, USA 

13 Department of Atmospheric, Oceanic, and Space Sciences, University of Michigan, Ann Arbor, Michigan, USA 

14 Department of Physics, University of Oxford, Oxford, UK 

'^Research Institute for Applied Mechanics, Kyushu University, Fukuoka, Japan 

16 ECMWF, Shinfield Park, RG2 9 AX, Reading, UK 

Correspondence to: B. H. Samset (b.h.samset@cicero.uio.no) 

Received: 22 October 2012 - Published in Atmos. Chem. Phys. Discuss.: 7 November 2012 
Revised: 8 February 2013 - Accepted: 12 February 2013 - Published: 1 March 2013 


Abstract. The impact of black carbon (BC) aerosols on the 
global radiation balance is not well constrained. Here twelve 
global aerosol models are used to show that at least 20 % of 
the present uncertainty in modeled BC direct radiative forc- 
ing (RF) is due to diversity in the simulated vertical profile 
of BC mass. Results are from phases 1 and 2 of the global 
aerosol model intercomparison project (AeroCom). Addi- 
tionally, a significant fraction of the variability is shown to 
come from high altitudes, as, globally, more than 40 % of the 
total BC RF is exerted above 5 km. BC emission regions and 
areas with transported BC are found to have differing char- 
acteristics. These insights into the importance of the vertical 
profile of BC lead us to suggest that observational studies are 
needed to better characterize the global distribution of BC, 
including in the upper troposphere. 


1 Introduction 

Unlike most atmospheric aerosols, black carbon (BC) ab- 
sorbs solar radiation. This warming effect of BC has led 
to suggestions, both in the scientific community (Editorial 
Nature, 2009; Grieshop et al., 2009; Shindell et ah, 2012; 
Hansen et ah, 2000) and among policy makers, for reduc- 
tion of BC emissions to mitigate global warming. The un- 
certainties in the radiative forcing (RF) of the direct aerosol 
effect of BC are however large (Ramanathan and Carmichael, 
2008; Schulz et ah, 2006; Feichter and Stier, 2012), hamper- 
ing mitigation studies (Koch et ah, 2011a). Among the causes 
of these model uncertainties are assumptions about the verti- 
cal concentration profiles of BC (Zarzycki and Bond, 2010). 
There are also significant discrepancies between models and 
observations (Koch et ah, 2009). 
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Aerosol radiative forcing is a measure of the effect a given 
change in aerosol concentrations has on the atmospheric en- 
ergy balance. The efficiency with which BC can induce RF 
is however dependent on external factors such as surface 
albedo, water vapor, background aerosol distributions and, 
most notably, the presence of clouds (Flaywood and Shine, 
1997; Zarzycki and Bond, 2010). The BC forcing is therefore 
highly sensitive to the full 3-D distribution of BC concen- 
tration (Samset and Myhre, 2011). Anthropogenic direct BC 
forcing estimated by the current major global aerosol mod- 
els varies between 0.05 Wm"’ and 0.38 Wm -2 (Myhre et 
al., 2013). This range is similar to earlier estimates (Schulz 
et al., 2006; Feichter and Stier, 2012), and difficulties in re- 
ducing the range highlights the urgent need to understand the 
different causes of this model variability in BC forcing. 

Several studies have previously indicated that the vertical 
transport is one area where the models still differ signifi- 
cantly (Koffi et al., 2012; Schwarz et al., 2010; Textor et al., 
2006, 2007). This study uses input from 12 global aerosol 
models participating in the AeroCom model intercompari- 
son project to compare and study the impacts of modeled BC 
vertical forcing profiles. By combining the models’ own con- 
centration profiles with a common 4-D (spatial and temporal) 
efficiency profile (EP) of RF per gram of BC, we here recal- 
culate and compare the exerted RF of the BC direct aerosol 
effect at various altitudes and spatial regions. By comparing 
calculations using the full 4-D in clear and cloudy skies with 
analyses using global and annual mean profiles, we can also 
isolate the contributions from the cloud field and variations 
due to regional differences. 


2 Data 

The AeroCom initiative asks models to simulate the direct 
aerosol radiative effect under as similar conditions as possi- 
ble. All models run single year simulations using the same 
base years for aerosol emissions (2000 and 1850), with un- 
changed meteorology to exclude indirect effects. Both circu- 
lation models and transport models have participated. Two 
sets of model intercomparisons have been performed, here 
labeled AeroCom Phase 1 (Dentener et al., 2006; Textor et 
al., 2006; Schulz et al., 2006) and AeroCom Phase 2 (Myhre 
et al., 2013). See the references or the main AeroCom web- 
site (aerocom.met.no) for details. 

For the present study, 4-D (spatial and temporal) BC con- 
centration fields from 12 global aerosol models participat- 
ing in AeroCom Phase 2 are used. For six of the models, 
concentration fields of BC from fossil and biofuel burning 
only (BCFF) are also used. The 1850 concentrations are 
subtracted from the 2000 ones, leaving the contribution due 
to anthropogenic BC or BCFF emissions. Unless otherwise 
stated, all results below are for anthropogenic aerosols only. 

Three of the same modeling groups also provided BC con- 
centration profiles to AeroCom Phase 1, and these fields 


are also included here. The IMPACT model has not under- 
gone major changes between the AeroCom phases, and for 
CAM4-Oslo the changes related to the treatment of BC are 
small. The OsloCTM2 model has however been heavily re- 
vised with the addition of both an ageing scheme for BC and 
improved treatment of BC washout between PI and P2. 

From AeroCom P2, participating models are NCAR- 
CAM3.5 (Lamarque et al., 2012), CAM4-Oslo (Kirkevag et 
al., 2013), CAM5.1 (Liu et al., 2012), GISS-modelE (Koch 
et al., 201 lb), GMI (Bian et al., 2009), GOCART-v4 (Chin et 
al., 2009), HadGEM2 (Bellouin et al., 2011), IMPACT (Lin 
et al., 2012), INCA (Szopa et al., 2012), ECHAM5-HAM 
(Zhang et al., 2012), OsloCTM2 (Skeie et al., 2011) and 
SPRINTARS (Takemura et al., 2005). From AeroCom PI, 
participating models are CCM-Oslo (Kirkevag and Iversen, 
2002) (earlier version of CAM4-Oslo, denoted UiO.GCM in 
AeroCom PI), IMPACT (Liu et al., 2005) and OsloCTM2 
(denoted as UiO CTM in AeroCom PI) (Myhre et al., 2003). 


3 Methods 

3.1 4-dimensional forcing efficiency profiles 

Samset and Myhre (2011) presented a set of vertical profiles 
of the aerosol direct radiative specific forcing (forcing per 
unit mass) for BCFF, i.e. the amount of top-of-atmosphere 
shortwave radiative forcing exerted per gram of aerosols at 
a given altitude. The profiles are for a particular model, the 
OsloCTM2, running under AeroCom P2 configuration. Here 
we term these efficiency profiles (EP). In the literature, the 
term normalized radiative forcing is also commonly used. 
We here employ the full 4-D efficiency profiles, either for 
all sky or clear sky conditions. For each model grid point and 
time step, we multiply the modeled BC concentration by the 
OsloCTM2 EP to get the contribution to the total shortwave, 
top-of-atmosphere BC radiative forcing. This gives us inter- 
comparable 4-D RF fields, something that is not immediately 
available from each model. 

Differences in model specific treatment of clouds, water 
uptake and microphysics are not treated by this method, i.e. 
they correspond to the ones for OsloCTM2. We are left solely 
with variations due to the concentration profiles and the total 
aerosol burden of the host models. We label this the recalcu- 
lated RF, to emphasize that it is heavily correlated with the 
burden and should not be taken directly as an estimate of BC 
forcing of each model. Since the model that was used to pro- 
duce the profiles of RF per gram has among the strongest 
global mean forcing efficiencies (Myhre et al., 2013), most 
recalculated RFs can be expected to be stronger than their 
host model would predict. However, by dividing the recal- 
culated RF by the total aerosol burden of the host model, 
we extract the variability in specific forcing that is due only 
to the vertical profile of BC produced by the host model in 
conjunction with the OsloCTM2 efficiency profile. We can 
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Fig. 1 . Modeled BC global mean (a) burden, (b) RF and (c) forcing efficiency. Yellow boxes indicate mean, one standard deviation and 
max/min values. Mean values and spreads for AeroCom PI and P2 (hatched whisker boxes) are taken from Schulz et al. (2006) and Myhre 
et al. (2013) respectively. 



subsequently use the vertical RF profiles to study this vari- 
ability (see Sect. 4.1). 

The remaining variability in specific forcing contains 
information on regional differences (high or low surface 
albedo, indigenous emissions, transport region or both, cloud 
fields), and aerosol optical properties and mixing tides. To 
separate these effects, we ran parallel analyses using only a 
global mean profile of RF per unit mass, a distinct profile of 
RF per unit mass for clear sky conditions, or a combination 
of clear sky conditions and a global profile. 

3.2 Region definitions 

A set of illustrative regions is chosen for some of the results 
and discussion below. “Europe” is defined as the box cover- 
ing longitudes —10 to 30, latitudes 38 to 60. “China” covers 
longitudes 105 to 135, latitudes 15 to 45. “Arctic” covers lati- 
tudes 70 to 90, all longitudes. Europe and China represent re- 
gions with high industrial BC emissions and significant frac- 
tions of the global BCFF emissions, making their total BC 
forcing sensitive to pure vertical transport and wet scaveng- 
ing. The Arctic is a region with low indigenous emissions, 
but important BC contributions at high altitudes transported 
from other regions. Its total BC forcing is therefore sensitive 
to model differences in vertical and long range BC transport 
and removal. 


4 Results 

4.1 RF variability due to vertical profiles 

We first quantify the fraction of the total modeled RF vari- 
ability attributable to vertical profiles alone. Figure la shows 
the global mean anthropogenic BC burden from models par- 
ticipating in AeroCom Phase 2, in addition to three models 
from AeroCom Phase 1 included to investigate the magni- 
tude of variability due to model developments. We find a 
multi-model mean of 0.19mgm -2 , but with a relative stan- 
dard deviation (RSD) of 32 % and a model spread from 0.09 
to 0.37 mg m -2 . Table 1 lists the numbers for individual mod- 
els. Under an assumption of equal specific forcing, this alone 
will cause a large diversity in the total BC forcing predicted 
by the models. AeroCom Phases 1 and 2 found a similar 
burden range, as shown by the whisker boxes in Fig. la. 
Note that the AeroCom PI includes BC from all sources, 
i.e. including biomass burning, while P2 results are for BC 
from fossil and biofuel burning (BCFF) sources only. This 
causes the P2 burden and RF means to be lower than for 
PL Figure lb shows the RF for each model, recalculated 
from the concentration profiles using the common EP from 
OsloCTM2 (Samset and Myhre, 2011). We find a stronger 
mean RF than from either of the AeroCom phases. This is 
partly due to the strength of the EP used, which is higher 
than for most AeroCom models, and partly to the fact that 
the present study includes BC from biomass burning as men- 
tioned above. See Discussion for further comments. 

The recalculated RF values are highly correlated with the 
burden values(Pearson corr. coeff. p — 95 %), as expected 
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Fig. 2. Comparison of modeled concentration and RF profiles, (a-c) BC concentration vertical profiles, global mean and for two selected 
regions. Overlain is the annual mean forcing efficiency profile for the selected region (grey dashed line). Solid lines show AeroCom P2 
submissions, dashed lines show PI. (d-1) BC RF per height, divided by the modeled global mean BC burden, globally and for three selected 
regions, (g-i) Vertical profile of integrated absolute BC RF. Lines indicate the 50 % mark and 500 hPa altitude. 


since a common EP is used. However, dividing the recalcu- 
lated RF by the global mean burdens gives global estimates 
of BC forcing efficiency (radiative forcing exerted per gram 
of BC aerosol), that are independent of the burden value sim- 
ulated by the host model. This is shown in Fig. lc. If the 
modeled spatial and temporal aerosol distributions, notably 
the vertical profiles, were also identical or merely related by 
a simple scaling, this spread would vanish. The remaining di- 
versity therefore carries information on the impact of the 4-D 
BC concentration profiles on the model RF spread, separate 
from the variability due to burden differences. 

In the following paragraph we quantify the contribution of 
the vertical profile to total model spread. Our method is to 
calculate the relative standard deviation (RSD) of the resid- 
ual variability, after normalizing the recalculated RF by the 
model burdens. We then compare this to RSDs from the na- 
tive model RF estimates. First, note that there is a resid- 
ual variability in the forcing efficiency (Fig. lc) that is due 
to the variations in the 4-D (spatial and temporal) aerosol 
mass distributions, with a relative standard deviation (RSD) 
of 16%. This number is dominated by the vertical distri- 
butions, but will also have components from horizontal and 


temporal variability in the aerosol fields between models. To 
isolate the contribution from vertical distributions, we per- 
formed the analysis using global mean BC and efficiency 
profiles per month, averaging out any horizontal differences, 
and then using global, annual mean profiles, to also remove 
temporal variability. The resulting RSD from vertical profiles 
alone is found to be 13 %. 

Next, we compare this to variability on forcing per bur- 
den in AeroCom PI and P2, which in both cases is around 
40%. Our residual RSD from vertical profiles is therefore 
approximately 30 % of the variability values from both these 
studies. However, three models in P2 have mass extinction 
coefficients that deviate significantly from recommendations 
given in the literature (Bond et ah, 2006) and thus unduly 
influence the RSD. Removing contributions from these three 
models reduces the RSD on the forcing per burden to 32 %, 
subsequently increasing our estimate of the variability due to 
vertical profiles to 40 %. 

Finally, we need to relate this to the total model RF vari- 
ability. Assuming that RF variability is in some way a combi- 
nation of variability in burden and forcing efficiency, we first 
note that in PI and P2 these two factors have approximately 
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equal RSDs. If they were uncorrelated the vertical distribu- 
tion could therefore be said to contribute half of the total vari- 
ability on RF, which from our above calculation is 20 %. This 
is however not the case. We will show below that 20 % can 
instead be considered a lower limit. 

To investigate the correlation of variability in burden and 
forcing efficiency, we note that in both AeroCom PI and P2, 
the burdens and forcing efficiencies are weakly anticorre- 
lated. This is apparent from the fact that the RSD on the RF 
(shown in Fig. 1) is lower than it would be if the errors on 
burden and forcing per burden were uncorrelated. Quantify- 
ing this effect using numbers from AeroCom P2, we find a 
weak Pearson correlation coefficient of p = —0.46 (or —0.40 
if we again remove the three outlier models). 

In the present analysis, we can estimate the impact of the 
vertical distributions by the fraction of mass simulated above 
500 hPa, or approximately 5 km, shown in Table 1 as M5k. 
M5k is calculated by integrating the BC vertical profiles pre- 
sented in Fig. 2a for model layers above or below an average 
pressure of 500 hPa, and then taking the ratio of these values. 
We observe that M5k is strongly correlated with the recalcu- 
lated RF (p — 0.70) and forcing per burden (p = 0.97), as ex- 
pected due to the efficiency profile used for the present anal- 
ysis. However M5k is also weakly positively correlated with 
the total burden (p = 0.48). The vertical distribution variabil- 
ity therefore does not contribute to the observed anticorrela- 
tion between burden and forcing per burden in AeroComPl 
and P2. 

We can assume that forcing per burden is detennined by a 
combination of the vertical profile (positive correlation with 
burden) and a set of uncorrelated global variables such as 
BC optical properties (which must then have a combined 
negative correlation with burden). The variability on BC RF 
would be higher if the models had not compensated for high 
burden with a low forcing efficiency and vice versa. Since 
the vertical variability rather leads to a high efficiency for a 
high burden, its impact on the RF variability is likely stronger 
than the minimum estimate of 20 % above. The present anal- 
ysis does not allow for rigorous quantification, but it is un- 
likely to be larger than the 50 % of the variability on total BC 
RF caused by forcing per burden. Hence we have a range of 
20 % to 50 % of the variability of modeled BC RF caused by 
differences in vertical distributions. 

4.2 Comparison of modeled RF profiles 

Next we compare global and regional BC vertical profiles 
and EPs. As also shown in previous model comparisons 
(Schwarz et al., 2010; Textor et al., 2006, 2007), there are 
large differences in the modeled BC concentrations. Fig- 
ure 2a-c shows the annual mean vertical profiles of BC 
concentration from all models, globally and for two se- 
lected regions (Arctic and China). Globally, the concentra- 
tions quickly decrease with altitude while the forcing effi- 
ciency increases. The Arctic is a region with negligible in- 


digenous BC emissions, but where the effects of transported 
aerosols are relevant (Shindell et al., 2008). While all mod- 
els show a significant BC contribution at high altitudes the 
model spread in the Arctic is significant, ranging from vir- 
tually constant below 200 hPa to double-peaked structures 
with maxima at 900 hPa and 200 hPa. This highlights the dif- 
ferences in transport and wet removal schemes between the 
climate models. AeroCom PI and P2 models can be seen to 
perform similarly where we are able to compare. China has 
large surface emissions of BC, and also sees contributions 
from transport at high altitudes, as is evident from Fig. 2c. 

Figure 2d-f shows the corresponding recalculated forcing 
as a function of altitude, divided by model layer height to 
remove differences due to the models’ vertical resolution. To 
highlight the vertical dependence, we have also normalized 
the profiles by the modeled global mean burden. The result- 
ing vertical profiles illustrate at what altitudes RF is exerted. 
Globally, RF can be seen to be mostly exerted in a range 
around 800 hPa for most models, with a secondary peak be- 
tween 400 hPa and 200 hPa for some models. The China re- 
gion follows the global pattern, but with the low altitude peak 
closer to the ground, reflecting the emissions there. The Arc- 
tic region, however, has a strong peak at high altitudes for 
most models, evidencing the relative importance of high al- 
titude BC RF there. These features are explainable by the 
differences in regional concentration and forcing efficiency 
profiles shown in Fig. 2a-c. 

Figure 2g-i illustrates the model spread in accumulated 
forcing, by integrating the absolute forcing per model layer 
from the surface and upwards. Globally, several of the mod- 
els exert 50 % of their forcing above 5 km; however there 
is a large spread around this value. Regionally, this picture 
is quite different. In the Arctic all models have a significant 
forcing component above 5 km, due both to a high fraction 
of aerosol at high altitudes and a strong forcing efficiency 
caused by high surface albedo. For the industrial regions of 
China the opposite is true, with the bulk of the aerosol located 
close to the ground where the forcing efficiency is weak. 

4.3 RF fraction exerted at high altitudes 

We find, in general, that a significant fraction of modeled BC 
RF is exerted at high altitudes, with a notable regional pat- 
tern. Figure 3 maps the model mean fraction of aerosol mass 
(Fig. 3a) and induced forcing (Fig. 3b) from model layers of 
altitude above 500 hPa, or approximately 5 km. Firstly, the 
distinction between BC emission regions and regions with 
mostly transported BC is evident. The former have small 
mass and RF fractions at high altitudes, typically 10-20%, 
while transport regions show RF fractions up to 80%. Fig- 
ure 3c compares the global and regional (Arctic, China and 
Europe) mean values of the mass and RF fractions. Note that 
the fraction of RF above 5 km is systematically higher than 
the mass fraction, due to the strongly increasing shape of the 
RF efficiency profiles. Globally, more than 40 % of the model 
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(a) Mass fraction above 5km 



180* I20*W 60*W 0* 60*E 120*E 180* 


Fig. 3. Black carbon mass and induced forcing at high altitudes, (a) Fraction of modeled BC mass above 5 km. (b) Fraction of modeled BC 
RF originating above 500 hPa (5 km), (c) Mean fraction of mass (orange) and RF (yellow) globally and for three selected regions. Boxes 
indicate one standard deviation on the model spread, whiskers show maximum and minimum values. 


simulated RF from BC comes from model layers above 5 km, 
while only 24 % of the mass is found in this region. This il- 
lustrates the importance of validating model vertical profiles 
and transport codes not only in industrial regions, but also in 
transport regions with low indigenous BC emissions. 

Figure 4 shows the RF and mass fractions in four al- 
titude bands, globally and for the three regions. Selected 
bands are up to 800 hPa, 800-500 hPa, 500-250 hPa and 
above 250 hPa, or approximately 0-2 km, 2-5 km, 5-10 km 
and above 10 km. Altitudes are chosen by finding the model 
layer with a global mean pressure that closest approximates 
the desired height. Such information is useful for comparing 
models with observational data, which are typically given for 
given altitude bands. This can in turn aid future constructions 
of best estimates for global BC forcing. We note that the pat- 
tern of high variability between model estimates of burden 
and forcing is apparent in all four bands chosen. 

4.4 Variability in forcing efficiency due to clouds and 
regional differences 

We have shown that, even on global mean, a significant frac- 
tion of the variability in the BC specific forcing is due to dif- 
ferences in vertical profiles when applying a common EP. It 
is instructive to further divide this variability into the compo- 
nents that make up the efficiency profiles. In Fig. 5a we show 
the zonal mean recalculated forcing efficiency for all models, 
and investigate the relative importance of the cloud field and 
of regional differences in albedo for the resulting spread. In 


Fig. 5c we have run the analysis using a global, annual mean 
efficiency profile instead of the full 4-D profile, and in addi- 
tion used clear sky conditions. The forcing efficiency has a 
weak but non-vanishing model spread, and is only weakly 
dependent on latitude since the effects of varying surface 
albedo, which normally strongly increases the forcing effi- 
ciency at the poles, are averaged out in the global profile. Not 
all vertical sensitivity of BC forcing is due to the aerosols 
being above or below clouds (Samset and Myhre, 2011). 
The remainder, caused predominantly by Rayleigh scatter- 
ing, water vapor and the competing effects of other aerosols, 
is the cause of the variability here. RSD on global mean val- 
ues is 10% (see Fig. 5b). Figure 5d shows the same analysis 
using the full 3-D efficiency profile. The RSD increases to 
13 %, and we now see the latitudinal effects of the high po- 
lar planetary albedo. Figure 4e shows the analysis using the 
global mean efficiency profile again, but under all-sky con- 
ditions. Again the RSD increases to 13%, due to the cloud 
field. The model variability in forcing efficiency that is due 
to vertical profile differences can therefore be decomposed 
into three factors: Equal and significant contributions from 
the cloud field and from regional differences, and a major 
contribution from the underlying sensitivity of BC forcing to 
altitude even in the absence of clouds and albedo differences. 
Harmonizing model treatment of clouds and albedo is there- 
fore not sufficient to remove uncertainties in BC forcing due 
to vertical profiles, highlighting the importance of ongoing 
aerosol model intercomparisons emphasizing other sources 
of variability (Stier, 2012). Combining the effects of clouds 
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Fig. 4. Mass and forcing fractions in four altitude bands, globally and for three regions. Upper rows show burden fractions in blues, lower 
rows show RF fractions in reds, (a) Global mean, (b) Arctic, (c) Europe, (d) China. 


and regional variations with the intrinsic vertical variability 
of the forcing efficiency yields the total variability picture 
seen in Fig. 5a. 

5 Discussion 

The model mean direct BC forcing recalculated above is 
stronger than in previous estimates (Schulz et al., 2006; Fe- 
ichter and Stier, 2012). This is mainly due to the fact that the 
host model for the efficiency profiles, OsloCTM2, has a high 
mean forcing efficiency compared to the others (Myhre et al., 


2013). The RF numbers in the present paper should therefore 
not be taken as new estimates for absolute direct RF forcing. 
We also observe that the recalculated RF values are highly 
correlated (95 %) with modeled burden values. This is again 
expected, and is due to the use of a common BC ER 

Different models will likely have different BC efficiency 
profiles, however we do not at present have enough infor- 
mation to estimate a model spread. A recent estimate using 
a column model (Zarzycki and Bond, 2010) shows forcing 
efficiencies above and below clouds comparable to the ones 
used in the present analysis. Flansen et al. (2005), using a 
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Fig. 5. Breakdown of spread in BC forcing efficiency, (a) Zonal mean of BC RF per unit mass using 3-D monthly resolved efficiency profiles 
and all sky conditions, for all models. Solid lines are results from AeroCom phase 2, dashed lines from phase 1. (c-e) Model mean (solid 
line) and maximum/minimum (dashed lines) when instead using (c) a global, annual mean efficiency profile and clear sky conditions, (d) 3-D 
profile and clear sky conditions, and (e) a global profile and all sky conditions, (b) shows the global mean values for the four cases. Boxes 
show one standard deviation on the model spread, whiskers show maximum/minimum values. 


global climate model, show (their Fig. 27a) a forcing profile 
that is of similar shape to the one used here, but with some- 
what weaker vertical gradient. Flowever while the absolute 
strength of the forcing efficiency will matter for the final RF 
estimates, only the shape, which is related to how clouds and 
microphysics are treated in individual models, will influence 
the variability due to vertical profiles alone. 

To attempt to quantify the sensitivity of the present analy- 
sis to the shape of the profile used, we reran the analysis with 
an EP that was weakened by 20 % at the top of the atmo- 
sphere (20 hPa) and unchanged at the surface, with a linear 
interpolation in between, resulting in an overall weaker EP 
gradient. This changed the global fraction of RF above 5 km 
by less than 5 %, indicating that the results are relatively sta- 
ble within reasonable variations of the EP. 

Results in the present study are given for total BC aerosol 
emissions only. Flowever, due to different source regions, the 
vertical profiles of BCFF could be different from BC. Six 
models also provided concentration profiles for BCFF. Us- 
ing the same efficiency profiles we performed the analysis 
also for BCFF, and found results consistent with what we 
have presented for BC (not shown). While the absolute forc- 
ing numbers differ due to lower total burdens for BCFF, the 
variability in vertical profiles is very similar to that for to- 
tal BC. Our conclusions here are therefore also applicable to 
model comparison results on BCFF only. 


Results presented here have all used emissions from year 
2000. One model (CAM4-Oslo) also provided simulations 
for year 2006. While the BC burden was 60 % higher for 
2006 emissions, the forcing efficiencies and vertical profiles 
were invariant. This gives confidence that the variability due 
to vertical profiles can indeed be regarded as independent of 
that due to the present day emissions dataset. 

Results for the PI and P2 submissions of the IMPACT 
model, where no major aerosol microphysical changes have 
been performed, are quite similar. For OsloCTM2, how- 
ever, where both an ageing scheme and modifications to the 
washout of BC were added, the forcing efficiency changes 
between PI and P2 by as much as half of the full range ob- 
served. Hence, the transport scheme and model treatment of 
BC are crucial factors in determining the modeled value of 
forcing efficiency, and both of these factors are closely linked 
to the vertical distribution. 


6 Conclusions 

We have shown that the previously documented large spread 
in BC aerosol concentration profiles is enhanced for vertical 
BC RF profiles. Using 12 global aerosol models, we show 
that most models globally exert 40 % of their BC forcing 
above 5 km, and that regionally the fraction can be above 
70%. The spread between models is however quite large, 
and we computed that at least 20%, possibly as much as 
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50 %, of the differences in modeled BC RF can be attributed 
to differences in vertical profiles. Harmonizing models with 
respect to the treatment of clouds, albedo and other factors 
relevant to the calculation of radiative forcing, is found to 
not be sufficient to remove this variability. To propose effi- 
cient mitigation measures for BC, its radiative forcing needs 
to be well understood both in emission and transport regions. 
Further model improvements and comparisons with data are 
needed, and should focus on both of these types of region. 
Observational studies are needed to provide input to model- 
ers, with the aim of better characterizing the global distribu- 
tion of BC, especially in the upper troposphere. It is however 
clear that while the BC vertical profiles are important, they 
are not sufficient to explain all the remaining differences be- 
tween global aerosol models. Further model intercomparison 
studies, e.g. on BC optical properties, surface albedo, treat- 
ment of clouds and aerosol transport and washout, are there- 
fore also needed if the impact of BC on the global radiation 
balance is to be sufficiently constrained. 
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